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ABSTRACT 

A direct observation of the spiral structure of the Galaxy is hindered by our position 
in the middle of the Galactic plane. We propose a method based on the analysis of the 
birthplaces of high Galactic latitude runaway stars to map the spiral arms and deter- 
mine their dynamics. As a proof of concept, the method is applied to a local sample 
of early-type stars and a sample of runaways stars to obtain estimates of the pattern 
speed (f2 P! iocai = 20.3 ± 0.5 kms _1 kpc _1 and f2 Pj runaway = 21.9 ± 8.6 kms _1 kpc _1 ) 
and the spiral arm's phase angle (0i OC al — — 53.1°±3.1° and Runaway = — 63.5°±22.5°). 
We also estimate the performance of this method once the data gathered by Gaia, in 
particular for runaway stars observed on the other side of the Galaxy, is available. 
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1 INTRODUCTION 

The spiral arms are readily visible in photographic images 
of face-on spiral galaxies as exemplified by the galaxy M51. 
They are usually traced, in visible light, by young hot stars 
and respective HII regions or, in other words, are the pre- 
ferred sites of star formation. Our position within the Galac- 
tic disc does not permit a direct observation of the spiral 
arms as they appear superimposed along our line of sight 
and are heavily obscured by dust. Nevertheless, many at- 
tempts have been made to determine their precise geometry 
and dynamics. The results obtained in different studies of 
the spiral structure of our Galaxy, using a variety of tracers 
(gas, dust an d stars) , were summarised and synthesised in a 
metastudy bv lValleel j 2008). The picture that emerges is that 
of a 4-arm spiral structure with pitch angle p — 12.8°. How- 
ever, some controversy remains regarding the spiral struc- 
ture properties, in particular the numb er of arms and their 
shape (pitch angle) : ILepine et al.l {2011), for example, claim 
the existence of "square" -shaped spiral arms, corresponding 
to the 4:1 resonance, instead of the usual logarithmic spiral 
arms. In any instance, only nearby segments of the spiral 
ar ms appear to be well de fined in most studie s (e.g. Fig. 3 
of Russcil 2003 and Fig. 4 of lLepine et al.ll201ll ). whereas the 
spiral structure behind the Galactic centre is virtually un - 
known (Zona Galactica Incognita according to Vallce 2005). 

The determination of the spiral pattern speed Q p is a 
closely related problem that has been the subject of many 
different studies. The differ ent methods us ed to estimate the 
value of n p are reviewed bv lGerhardl (|201ll ). The most direct 
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method consists in finding the birthplaces of open clusters 
or individual young st ars by computing their orbits (e.g . 
Amaral & Lepine 1997; Fernandez. Figueras & Torra 2001: 



Dias fc Lep ine 200|; iNaoz fe Shavivll2007l ). 



In this paper we propose the use of a different tracer of 
the spiral arms structure: the Galactic population of high 
Galactic latitude runaway stars. It has already been shown 
that high quality astrometry, in particular that obtained 
with Hipparcos (see ?), permits the recovery of the birth- 
places of runaway stars thus, the higher quality promised 
by the Gaia mission, will allow the analysis of a much 
larger volume of the Galaxy. Although it is more common 
to use large r structures a s tracers, as the af orement ioned 
HII r egions jRusseilll2003l ) or open clusters |Dias fc Lepind 
2005), sin gle stars were also previously used (e.g. Chepeid 
stars by iMaiaess. Turner fc Lane] 120091 and O-B stars by 
iFernandez et al l 200ll ). However, "normal" early-type stars 
are confined to the disc and so are heavily affected by inter- 
stellar reddening, whereas high Galactic latitude runaway 
stars are not. Moreover, since runaway stars can travel large 
distances (several kpc) away from their birthplaces, it is pos- 
sible to map portions of the spiral arms that are further away 
than the runaway stars. 

Nevertheless, the number of high Galactic latitude run- 
away stars and the spatial distribution of their birthplaces 
does not permit yet the application of this method with 
high accuracy, however the number of known stars of this 
type will increase by at least one order of magnitude with 
the Gaia satellite mission, which will deliver high accuracy 
astrometric data. Thus, in this paper our objective is to con- 
duct a proof of concept study to show how the use of the 
birthplaces of runaway stars may be used to trace the spiral 
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arms so it can be applied to the larger population of runaway 
stars observed by Gaia. 

In order to accomplish this objective we apply a vari- 
ation of the method used in other studies (usually with 
open clusters as tra cers, e.g. lAmaral fc Lepinel 1 19971 and 
iDias fc Lepinel l2005h . This method consists in using kinc- 
matical information about an object to compute its orbit 
in the Galaxy's potential and thus its birthplace (given an 
age estimate). Since objects with different ages indicate the 
position of the spiral arms in different time instants, this 
method also provides an estimate of the pattern speed Q p . 
Because of the shortcomings of the runaway stars sample we 
apply this method to a sample of local early-type stars as it 
provides a more accurate test (since it is larger and we may 
use the more accurate Hipparcos astrometry). After show- 
ing the results obtained with these two samples (local and 
runaways), we discuss the expected improvement in perfor- 
mance when this method is applied to Gaia data. 



2 PATTERN SPEED DETERMINATION 
2.1 Samples 

The sample of runaway stars un der consideration in t his pa - 
per is the sample presented in I Silva fc Napiwotzkil (|201lh . 
This is a sample of 96 stars covering the brightness range 
6.5 < V < 14.5, and an altitude above the Galactic plane of 
0.3 — 30.5 kpc. The ejection velocities, fli ght times and birth- 
places were computed as detailed in ISilva fc Napiwotzkil 
i|201ll ). In summary, the orbits were integrated backwards in 
an axisymmetric gravitational potential. Thus, the influence 
of the slight overdensities in the region of the spiral arms 
(less than 10 per cent of the total gravitational potential, ac- 
cording to lFernandez et ai1l200ll ').and the Galactic bar were 
ignored. This is justified for the following reasons: the total 
effect of the spiral arms in the velocity field is estimated 
to be only about 3 kms -1 , according to iFernandez et al.l 
l|200ll ), which corresponds to a maximum systematic effect 
in the position of the birthplaces of ~ 600 pc for the older 
stars (but less for the younger ones); this method is meant 
to be used with high Galactic latitude runaway stars and it 
was already shown in lSilva fc Napiwotzkil (|20 1 if ) that these 
stars are not very sensitive to changes in the disc Galac- 
tic potential, given that they leave the Galactic plane with 
large ejection velocities. The 51 birthplaces which were de- 
termined with a precision < 1.5 kpc in both x and y direc- 
tions are shown in Fig. [1] 

In principle it would be possible to fit directly these 51 
points to determine the feasibility of the proposed method 
to determine the spiral arms shape, position and pattern 
speed, however this sample is small and the distances and 
proper motions have relatively large measurements errors 
(no parallaxes). For this reason a second sample, consisting 
of local early-type main sequence stars, was selected allowing 
for a better determination of birthplaces. This sample and 
the runaway stars sample will be referred to as sample A 
and sample B, respectively. 

Sample A was selected fr om a cross-match o f the 
Hipparcos (new reduction by van Leeuwenl 120071 ) and 
the Bright Star |Hoffleit fe Jaschekl Il99lf ) catalogues. The 
selection criteria was: B — V colour < 0.05 (corre- 
sponding to an effective temperature < 9500 K, cf. 
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Figure 1. Birthplaces of the runaway stars in the Galactic plane 
(frame centered in the centre of the Galaxy) including 1 a error 
bars. The concentric circles indicate distances of 5, 8 and 10 kpc 
to the centre. 




Figure 2. Deredenned colour-magnitude diagram of local early- 
type stars (sample A) with respective error bars. The ev olution- 
ary tr acks obtained through the bolometric corrections of lFlowei] 
l|1996l ) from the theoretical tracks by ? are also plotted. 



iNapiwotzki. Schonberner fc Wenskelll993f ). and parallax > 
5 mas (the Bright Star catalogue is limited to a distance of 
~ 200 pc for stars with Mv = 0) . After the removal of stars 
with spectral type later than B (according to the Bright Star 
catalogue classification) and of spectroscopic binaries, the 
sample was reduced to 516 stars. Note that the colours were 
deredenned using a procedure base d on Stromgren uvby/3 
photometry l|Napiwotzki et aL 19931). The uvbyf] pho tome- 
try was obtained from the Hauck fc Mermilliod (|l998T ) cata- 
logue. In Fig.[2]the colour-magnitude of sample A is shown. 



2.2 Age and velocity determination 

The age of each individual star in sample A was esti- 
mated using a method analogous to the one used to com- 
pute the evolutionary ages of the runa way stars (sample B) 
ijSilva fc Napiwotzkil l201ll . also see eg. Ilrrgang et al.l l201fj| : 



iRamspeck et al]|200ll ). However, here the most recent evolu- 
tionary tracks by ? for non-rotating stars were used instead 
of the older models by the Geneva group. These tracks were 
computed using updated input physics and the revised value 
of the solar metallicity (?). The tracks were first converted 
from the L — T e g to the My — (B — V)o space and then 
linearly interpolated in order to obtain an evolutionary age 
estimate. 

The conversion of the theoretical tracks was accom- 
plished in two steps: the absolute magnitudes in V band 
corresponding to the given lu minosities were computed using 
the bolometric corrections by Flower (1996); then the (B — 
V)o colours corresponding to the given effective tempera- 
tures were computed using the formulae in lNapiwotzki et al.l 
(|l993h - These converted tracks are shown in Fig. [21 

The full space velocities were computed for all stars 
in both sample A and sample B were computed from the 
distances, measured proper motions and radial velocities. In 
the case of sample A, the distances were derived from the 
respective Hipparcos parallaxes. 

2.3 Fitting procedure 

The method adopted to compute the pattern speed consisted 
in two steps: 1) tracing back the orbit of each star in order to 
determine its birthplace (assumed to lay on a spiral arm); 
2) rotation of the position of the birthplace by the angle 
corresponding to the pattern speed fi p and the age of the 
star. The best estimate of the pattern speed is the one that 
minimises the distance of the birthplaces to a given model 
of the spiral arms structure, for the sample. Note that other 
parameters of the model may be determined if a suitable 
sample is provided. Thi s is e ssentially the method used by 
I Amaral fc Lepind l| 1997ft and lDias fc Lepind l|2005l ) on their 
sample of open cluster s. In our cas e we have adopted the 
cartographic model by IValled | 20081 ) as a representation of 
the present spiral arms structure as it based on previous 
studies using a variety of tracers and it has a simple math- 
ematical representation (logarithmic spiral). Thus, in polar 
coordinates: 

r(0) = r e bB+nn/2+4 ', (1) 

where b — tanp, <f> is the initial phase or phase angle, 
and n = 0, 1, 2,3. The pa rameters of the spiral structure, 
as derived by[Vallcc (2008]), are: the pitch angle p — 12.8°, 
the initial radius ro = 2.1 kpc, and the initial phase <f> = 
—20°. However, these parameters were obtained assuming a 
distance of the Sun to the Galactic Centre, Rq — 7.6 kpc, 
whereas in our study we assumed Rq = 8 kpc. Taking this 
difference into account, a fit of the spiral s tructure to the 
tangent directions indicated in Vallee ( 2008) gives an initial 
phase <f> = —40° ■ • ■ — 45°, for Rq = 8 kpc keeping the other 
parameters fixed. Given the nature of our samples (sample A 
only includes local stars and sample B is small), from these 
parameters only the initial phase <f> was fitted together with 
the pattern speed Q p . 

Thus, we used the following target function in our min- 
imisation procedure: 

^^) = E^^^w, (2) 

where 
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dn(i) = s] M^p, u) - x B (e n , 0)) 2 + G/i(n p , to - y s (e n , </>)) 2 (3) 

is the minimum distance between the birthplace ro- 
tated by Q p (with coordinates Xi,yi) and the closest spiral 
arm (with the closest point having coordinates x s ,y s corre- 
sponding to an angle 9 n ), t\ is the estimated age of the star, 
and <7x(i) and o y (i) are the uncertainties in the coordinates 
of the birthplaces as derived through a Monte Carlo error 
propagation procedure from the uncertainties in the input 
observable parameters. Since the fitting procedure assumes 
the spiral arms have no thickness (they are treated as lines) 
the uncertainty in the position of the stars was assumed 
to be a minimum of 200 pc (corresponding to a 3 j thick- 
ness of 1.2 kpc, cf. iLepine et al.ll201ll : rRusseilll2003l ). so the 
birthplaces far from the centre of the spiral arms are not 
excessively penalised by small weights during the minimi- 
sation procedure. This also alleviates the problem of stars 
born outside the spiral arms having a large influence in the 
fitting with out resorting to some for m of sigma-clipping a s 
was done bv lNaoz fc Shavivl (|2007h or lDias fe Lepind l|2005l) . 
for example. 

Since the target function, _F(f2 p , </>), does not strictly 
follow a x 2 distribution, and hence does not have known 
statistical properties, the error in the derived parameters 
was determined through a Monte Carlo procedure. Thus, 
2000 replicas of the original sample were created, varying 
the input parameters (colour, parallax, proper motion and 
brightness) according to the respective errors distributions 
(assumed Gaussian) . As new age estimates are computed in 
each realisation of the Monte Carlo procedure, this technique 
is also a good way of propagating the age uncertainties. 



3 RESULTS AND DISCUSSION 
3.1 Pattern speed 

Although in principle it would be possible to constrain all 
the parameters of the simple logarithmic spiral model given 
by Equation fT}, in the concrete case of our samples this is 
not possible as the tracers used either cover only a very small 
part of the Galaxy (sample A), or are too few with large 
measurement errors and small spread in ages (sample B). 
Thus, only a fit to the phase angle, together with the pattern 
speed, was attempted. 

In order to look for the minimum of the target func- 
tion (Equation [2} a grid spanning the area of interest of the 
parameter space was prepared in order to find the location 
of possible solutions which were later more precisely deter- 
mined using an implementation of the Nelder-Mead sim- 
plex algorithm in pythoisQ. In Fig. [3]the contour plot corre- 
sponding to sample-A is shown. Only one solution appears 
to exist (within the equivalent to 6 a) which is repeated for 
phases separated by 90°, i.e. when one arm switches to the 
next one. However, care should be taken with contamination 
by stars formed in regions un related to the spir al arms, e.g. 
the Gould Belt. According to lTorra et~aH |2000l) . the Gould 
Belt has an age of ~ 60 Myr, thus in order to remove this 
feature we computed a new grid considering only stars older 

1 http:/ /docs. scipy.org/doc/scipy/reference/generated/ 
scipy.optimize.fmin.html 
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Figure 3. Contour plots of the x 2 values for sample A stars. 




Figure 4. Contour plots of the \ 2 values for sample A stars older 
than 60 Myr (to partially remove the Gould Belt). 

than ~ 60 Myr. The resulting contour plot is shown in Fig. 0] 
and it is visible that although there is still only one solution 
it appears to become more uncertain. 

The reason for this increased uncertainty is understood 
when the same procedure is applied to a sample of even 
older stars, as shown in Fig. [S] In this case the solution be- 
comes degenerate, showing that a sample of older stars does 
not provide a good constraint on the phase angle. On the 
other hand, the sample of stars younger than 60 Myr, shown 
in Fig. [6] provides a solution that behaves in the opposite 
way, providing a better constraint on the phase angle but 
performing worse in constraining the pattern speed. This is 
simply a reflection of the obvious fact that younger objects 
are good at describing the present state of the spiral arms 
but have not travelled enough to be able to carry informa- 
tion about the pattern speed, whereas the opposite is true 
for older objects. 

In fact, if the same procedure i s appl ied to the young 
clusters (< 7 Myr) of the lPias et alj <|2002l > sample, we find 
that the present phase angle appears closer to —50° than 




Figure 5. Contour plots of the \ 2 values for sample A stars older 
than 90 Myr. 




Figure 6. Contour plots of the \ 2 values for sample A stars 
younger than 60 Myr. 

to the -20° derived by I Vailed (|2008i ) . However, as we have 
seen, most of this difference is attributable to a difference 
in the chosen distance of the Sun to the Galactic Centre 
(Rq = 8 kpc in our case), and, in particular, the constrains 
used bv lValleel (|2008ft imply ~ -45° assuming Rq = 8 kpc. 
In Fig. [7] we have plotted the grid of the parameter space 
for this sample of clusters, and in Fig. [8] a comparison be- 
tween the two different angles and how they compare with 
the sample of young clusters. Although a phase angle of 
—20° appears to be a better fit in some regions, a larger 
angle fits better other regions (Perseus arm, for example). 
Even though th e agreement b etween our solution and the 
one obtained by IValleel (|2008h (and hence with constraints 
provided by the tangential directions and t he distance to 
the Perseus arm obtained bv lXu et alj|2006l) is not perfect 
it would be possible to find a solution consistent with the 
constraint, however that would require further analysis be- 
yond the scope of this paper. Note that in the plot shown in 
Fig. [7] a phase angle of —20° is actually in the region of the 
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Figure 7. Phase angle estimate fr om sample of clusters younger 
than 7 Myr Cfrom lDias et alj|2002h . 




Figure 8. The spiral arms model bv lValled l|2008h with = -20° 
(solid line) and </> = —45° (dashed line), derived for Rq = 7.6 kpc 
and Rq = 8 kpc respectively, and <j> = —55.5° (dotted line). The 
different models are compared with the sample of open clusters 
younger than < 7 Myr (blue dots). 



worst solutions, but a phase angle of —45° (the corrected 
value for Rq = 8 kpc) is in the region of best solutions. 

The minimisation of the target function (Equation [2]) 
around the solution found in the parameter space, (fl p = 
20, <f> = —50°), yielded the values seen in Table [T] for three 
different samples: all of sample A, only stars in sample A 
that are older than 60 Myr, and the sample B. Note that 
the solution in the case of sample B is quite degenerate, 
however there is a solution compatible with the one found 
in the other two cases. 

The spiral arm structure with the phase estimated from 
sample A is shown in Fig.[9l together with the birthplaces of 
stars from samples A and B rotated from their original po- 



Table 1. Pattern speed (f2 p ) and phase angle (<j>) obtained for 
three different samples: entire sample A (local Hipparcos stars), 
only stars older than 60 Myr in sample A, sample B (runaway 
stars). 



Sample 


Q p (kms 1 kpc ) 


<t>n 


A (all) 


20.3 ±0.5 


-53.1 ±3.1 


A (> 60 Myr) 


21.2 ±1.1 


-60.3 ±7.4 


B 


21.9 ±8.6 


-63.5 ±22.5 




Figure 9. Distribution of the birthplaces in their present posi- 
tions, resulting from the derived pattern speed. Sample A and 
sample B (runaway stars) are plotted as coloured circles (with 
sizes proportional to the uncertainty in the position) and squares 
with error bars, respectively. Only runaway stars with standard 
deviations < 1.5 kpc are shown. The spiral arms are also plotted 
in their current position 



sitions by the derived (from sample A) pattern speed times 
the ages of the stars. It is interesting to note how the run- 
aways seem to have been born in many different arms, and 
in particular how well they appear to follow the "yellow" 
arm (to be identified with the Perseus arm). 



3.2 Comparison with other studies 

Estimates of the pattern speed in the literature cover a 
wide range of values: using different variations of the "birth- 
place" determinat i on me thod with samples of open clusters, 
Amaral fc Lepine (Il997l ) found f2 p ~ 20 kms -1 kpc -1 and 
Dias fc Lepind l|2005l ) O r « 24±1 kms " 1 kpc" 1 ; by fit- 
ting a kinematical model of the Milky Wav lFernandez et al.l 



l|2001f) found S7 P w 30 kms" 1 kpc" 1 (with an uncertainty 
between 2 kms" 1 kpc" 1 and about 7 kms" 1 kpc" 1 depend- 
ing on the specific tracer); iMartos et al.l (|2004T ) obtained 
Q p ~ 20 kms" 1 kpc" 1 by fitting dust obse rvations with a 
dynam ical model of the Galaxy. Moreover, iNaoz fc Shavivl 
(2007) estimate different pattern speeds, using a variation 
of the "birthplace" method with a sample of open clus- 
ters, for different structures (with an uncertainty close to 
~ 2 kms" 1 kpc" 1 for all estimates): the Sagittarius-Carina 
arm is actually a superposition of two spiral structures with 
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Pl i ~ 16.5 kms _1 kpc _1 and £l Pj 2 ~ 29.8 kms -1 kpc -1 , 
the Perseus arm with f2 p w 20 kms -1 kpc -1 , and the 
Orion "armlet" with Q p ~ 28.9 kms - kpc -1 . By fitting 
models of gas flow i n the Milky Way to CO observations, 
iBissantz etafl l|2003h obtained fi p w 20 kms" 1 kpc -1 . 

The different estimates of the pattern speed seem to 
suggest (even though it is not clear how uncertain some 
of these estimates are) a value in the r ange fi p ~ 20 — 
25 kms" 1 kpc- 1 (see also iGerhardl 1201 if ), which would be 
compatible with our own estimates. Furthermore, there is 
independent confirmation, from the observation of the pre- 
dicte d ring-shaped gap i n the radial HI density distribu- 
tion l|Amores et al.ll2009f ). of a corotation radius close to 
the value of the distance of the Sun to the Galactic cen- 
tre, implying a pattern speed Q p ~ 25 — 29 kms -1 kpc -1 . 
As we have seen, our best estimate gives a value of f2 p = 
20.3 ± 0.5 kms -1 kpc -1 for the pattern speed, which, al- 
though consistent with other determinations, appears to be 
at odds with estimates of the corotation radius, particularly 
given the small uncertainty. However, we must bear in mind 
certain caveats to this result: the uncertainty includes only 
the formal statistical contributions, ignoring systematic ef- 
fects; it is possible that the inclusion of regions of star for- 
mation unrelated with the spiral arms may introduce biases 
in the estimates. Indeed, it can be seen in Table [T] that re- 
moving most of the Gould Belt from the sample produces 
an estimate (fi p = 21.2 ± 1.1 kms -1 kpc" 1 ) closer to the 
upper limit of the suggested range for the pattern speed 
(f2 p ~ 20 — 25 kms -1 kpc -1 ), albeit with a larger uncer- 
tainty caused by the loss of the young stars constraining the 
current position of the arms. 

For the same reason, it is possible that the uncertainty 
on the phase angle determination was underestimated, be- 
cause the large number of young stars in the solar neigh- 
bourhood may introduce a bias in the position of the spiral 
arms, as they might have been formed in other regions (like 
the Gould Belt, or interarm features). However, the values 
we obtained for the phas e angle are close to the estimates 
bv lFernandez et all (|200lf ) (varying between cj> = -45° ±32° 
and cj> = 8° ± 52°, when assuming models with four spiral 
arms), and also with the estimate obtained from the sample 
of open clusters (see Fig. [SJ). 



4 PROSPECTS WITH GAIA 

The Gaia mission will deliver trigonometric parallaxes, im- 
proved proper motions and spectral energy distributions for 
many early-type runaway stars. This will: 1) improve the 
achievable accuracy for the existing sample, and 2) allow 
the selection of a much increased sample of runaway stars. 

In order to illustrate the expected increase in perfor- 
mance of our method once Gaia data becomes available, we 
have determined the birthplaces of two representative stars 
of our sample of runaways, one far away and one nearby, 
assuming the knowledge of their parallaxes and proper mo- 
tions with the predicted Gaia accuracy. The results obtained 
are summarised in Tableland contrasted with the previous 
determination. Our results for the nearby star HIP 59955 are 
already fairly accurate. Gaia measurements will reduce the 
overall uncertainty of the determination of the birthplace to 
a level comparable with the typical width of a spiral arm. On 



Table 2. Precision of the determination of the birthplaces of 
the stars EC 04420-1908 and HIP 59955 using present data 
l|Silva k. Napiwotzkil201ll) and the precision predicted using Gaia 
data. Estimates of the distances and the coordinates of the birth- 
places obtained using present data are indicated with a p sub- 
script whereas the estimates obtained using Gaia precision data 
are indicated with a g subscript. 







EC 04420-1908 


HIP 59955 


V 




13 


9.71 


dp 
J g 


(kpc) 
(kpc) 


io.44±*;»g 

10.55 ± 1.2 


i cr+0.38 

^•"-o.so 
1.5 ±0.02 


*p 
*g 


(kpc) 
(kpc) 


-11.4311; 1 ? 
-12.14 ±0.88 


. 31 +l.a7 
-4.30 ±0.60 


>p 


(kpc) 
(kpc) 


-7 40+ 2 - 74 
-8.22 ±0.78 


-4.39 ± 0.34 



the other hand, the improvement for the star EC 04420-1908 
is dramatic, with the X coordinate uncertainty decreasing 
by ~ 500 per cent and the Y coordinate uncertainty decreas- 
ing by ~ 300 per cent. Although the uncertainty obtained in 
this case is larger than the typical width of a spiral arm, it 
is no w smaller than the expected interarm distance 

Noting that the distance uncertainty that would be ob- 
tained from the Gaia parallax in the case of EC 04420-1908 
is ~ 1 kpc, it is clear that this error should be considered 
an upper limit if we hope to obtain precise estimates of 
the birthplaces' coordinates. The parallax errors that cor- 
respond to an error of 1 kpc in the distance, as a function 
of distance, are plotted in Fig. 1101 together with the paral- 
lax errors for stars with V < 13 and V < 14 (according to 
the formula in the Gaia Science Performance websit^El)- It 
is possible to conclude from this plot that it will be possible 
to determine directly - with high precision - the birthplaces 
of stars with V < 13 and V < 14 as far as 10 kpc and 
8 kpc, respectively. For stars further away other methods 
will be needed in order to attain a high precision, however 
it is important to note that distances of 14 kpc and 12 kpc, 
respectively for stars brighter than V = 13 and V = 14, will 
provide parallax errors of « 15 per cent. 

4.1 Sample selection and contamination with old 
evolved stars 

An issue always present in studies of runaway stars is the 
contamination by low mass stars in advanced stages of evo- 
lution, in particular the post-AGB and core Helium burning 
Horizontal Branch, as they may mimic the a t mospheric pa - 
rameters of main sequence stars (|Tobinll 19871 : lMartirJl2004h . 
However, as a direct measurement of the distance will be 
available with Gaia it will be possible to select runaway 
stars based on their absolute magnitudes My and effective 
temperatures T e g. 

From Fig. 1101 we have already seen that it will possible 

2 http://www.rssd. esa.int/index.php?project=GAIA& 
page=Science_Performance 
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> 0.15 




Figure 10. Parallax error as a function of distance corresponding 
to an error of 1 kpc (solid line), 2 kpc (dashed line) and 3 kpc 
(dotted line) in distance. The shaded areas correspond to the 
errors for stars with G < 15 (dark grey), G < 14 (medium grey) 
and G < 13 (light grey), that are predicted for Gaia. 



Figure 11. Comparison between the expected relative error, as 
a function of distance, for the distance estimates: from Gaia par- 
allaxes for a G < 13 star (solid line), from spectroscopy when 
the uncertainty in l ogg i s 0.05 dex (dashed line), as obtained by 
iNieva fc Przvbillal ||2012| 1. taking all error sources into account 
(dotted line). 



to determine the distance with high precision for stars with 
G < 13, and, more generally, with an error < 15 per cent 
for stars with G < 14, at least for distances up to 12 kpc. 
However, it can also be seen that an error < 30 per cent will 
be possible for stars with G < 15. Thus, considering that 
this error in the distance translates to an error of 0.65 in 
the absolute brightness My, it will be possible to use Gaia 
parallaxes to select a sample of early type runaway stars, 
even in those cases where the precision is not good enough to 
be used directly, as can be seen in Fig. [T^] This is because the 
absolute magnitudes of stars with spectral type B7 or earlier 
(effective temperature ~ 13500 K) are much lower than the 
absolute magnitude of the Terminal Age Horizontal Branch 
(TAHB, the end of Helium burning in the core), which is 
Mv > 1, for the same temperature. Note that this allows 
the selection of B7 stars as far as 12 kpc. 

Once the sample is selected, better estimates of the dis- 
tance could be obtained after a careful spectroscopic analysis 
in the case of the distant stars. Since the limiting factor of 
spectroscopic distance estimates is usually the uncertainty in 
the log g estimate, it is instructive to see how the error that 
is expected from Gaia parallaxes compares with good spec- 
troscopic es_timates 2 _with^maJluncertaint in log p. For 
example, INieva fe Przvbillal l|2012h obtain a distance accu- 
racy of ~ 10 per cent in their study of nearby B stars, from 
an analysis where the uncertainty in the log<? determina- 
tion is ~ 0.05 dex. The comparison between these estimates, 
and also the expected error if only the log g uncertainty (of 
~ 0.05 dex) is considered, is shown in Fig. 1111 Note how the 
Gaia estimates have acceptable errors only up to distances of 
10 to 12 kpc. Moreover, since the radial velocities measured 
with Gaia will have acceptable errors (< 10 kms -1 ) only 
in the case of the brightest stars (V = 13 . . . 14), follow-up 
spectroscopy may be needed in order to obtain good esti- 
mates of the radial velocities of the more distant stars. 

Also important, for both the selection of runaway stars 
and the follow-up spectroscopic analysis in the relevant 
cases, are good estimates of the effective temperatures. 



These could be obtained from Gaia spectrophotometry, if 
we assume it will be at least as accurate for effective tem- 
perature determination as medium band photometry, e.g. 
Stromgren- Crawford uvbyj3 phot ometry. The expe c ted er - 
ror would be then, according to iNapiwotzki et all l|l993h . 
about 3 per cent in the range 11000 K «C T eS 20000 K 
and 4 per cent for hotter stars. These should be conserva- 
tive estimates of the errors, which cou ld be further refine d 
by adding Galex UV photometric data dMartin et al.j|2005l ) . 

As was already suggested, the expected contamination 
by Horizontal Branch stars, for stars with spectral type B7 
or earlier will be extremely low, if the appropriate selec- 
tion criteria is adopted. Adopting a cr iteria similar to the 
one used by ISilva fc Napiwotzkil l|201ll ) (dependent on es- 
timates of logy obtained from spectroscopy) should limit 
the contamination to less than 10 per cent of the total. In 
fact, most contamination by Horizontal Branch stars can be 
eliminated by applying a cut in the T e g — Mv space (e.g. 
T cff > 13500 and M v < 0), as can be seen in Fig. [12] How- 
ever, this can be improved by combining good distance esti- 
mates with log g information, since the distances of 0.5 Mq 
and 2.5 M© stars of the same observed brightness differ by 
\/5. The contamination by post-AGB stars is more difficult 
to ascertain, however these stars are very rare and usually 
account for a very small amount of contamination in studies 
of ru naway stars (e.g. ~ 8 per cent in t he study by Martin 
20061 . and ~ 5 per cent in the study bv lSilva fc Napiwotzki 
201 II ). Nevertheless, the proposed criteria can be applied to 



these stars and also to stars in other post-Horizontal Branch 
evolutionary stages. 

In summary, combined with a small proportion of Hori- 
zontal Branch stars evolved past the TAHB, the level of con- 
tamination with old evolved stars should be < 10 per cent. 
Furthermore, we will be able to use Gaia astrometry to se- 
lect runaway stars with spectral types earlier than B7 up 
to a distance ~ 14 kpc and spectral types earlier than B3 
up to a distance ~ 22 kpc, assuming a brightness V = 15. 
However, follow-up spectroscopic analysis will be needed to 



8 M. D. V. Silva and R. Napiwotzki 




r ert (K) 

Figure 12. The Zero Age Horizontal Branch (ZAHB) and the 
Terminal Age Horizo ntal Branch (TAHB), from the models by 
iDorman et al.l Jl993) for stars with an Helium mass fraction 

Y = 0.247 and a metallicity [Fe/H] = —1.48, and the theoretical 
ev olutionary tra c k of a low mass Post-AGB star, from the models 
by ISchonbernen jl979h . in a My — T c g diagram. The absolute 
magnitudes for the Horizontal Branch models were o btained by 
applying the bolometric corrections of Flower (199^). The red 
dots corr espond to sta r s iden tified as Horizontal Branch in the 
survey by I Brown et al.l d200ST) . 

obtain good estimates of the distances in the case of the 
fainter stars. On the other hand, for stars brighter than 

V = 13 (corresponding to a distance ~ 9 kpc for B3 and 
~ 5 kpc for B7 stars) it will be possible to obtain good 
distance estimates directly from Gaia parallaxes. 

4.2 Spiral arms on the other side of the Galaxy 

As we have seen, it will be possible to select a sample of run- 
away stars reaching distances up to ~ 14 kpc, creating the 
opportunity to study the spiral arms behind the centre of 
our Galaxy. It is thus important to have an estimate of how 
many runaway stars will be observed behind the Galactic 
centre. The estimated number density o f runaway stars high 
abov e the Galactic plane is ~ 9 kpc -2 (jSilva fc Napiwotzkil 

2011) . This corresponds to a total of ~ 6400 stars in the 
Gala xy, in a 15 kpc radius (see also iNapiwotzki fc Silval 

2012) . Of these, about one third are at a distance > 3 kpc 
from the Galactic plane. 

From the reddening maps of ISchlegel et al.l (|l998l ) it is 
possible to estimate the expected extinction for stars out- 
side the Galactic plane, and in particular to obtain an upper 
limit to the expected extinction for runaway stars high above 
the Galactic plane, considering that the greatest amount 
of extinction occurs in the direction of the Galactic cen- 
tre (I = 0°). The extinction for stars 3 kpc away from the 
Galactic plane depends on the distance under consideration, 
as the Galactic latitude b decreases for larger distances. In 
Table [3] we show the values of extinction Ay obtained from 
the Schlegel maps for distances in the range 10 — 16 kpc. A 
distance of 10 kpc should be enough to probe the Sagittarius- 
Carina arm, 12 kpc should probe the Perseus arm, 14 kpc 
should probe the Cygnus arm, 16 kpc should probe the Crux- 



Table 3. Extinction predicted by the lSchleeel et ah red- 
dening maps for different Galactic latitudes. The Galactic lati- 
tudes were computed for points 3 kpc above the Galactic plane 
and at the given distances (d). Also indicated are the apparent 
magnitudes for stars with spectral types B3 and B7 when red- 
dened by the corresponding amounts. The values for the Northern 
and Southern Galactic hemispheres are distinguished by the use 
of parentheses, with the values inside parentheses corresponding 
to the Southern hemisphere. 



d (kpc) 


w n 


A v 


V (B3) 


V (B7) 


10 


17 


1.7 (0.4) 


15.6 (14.3) 


16.7 (15.4) 


12 


11 


1.0 (0.4) 


15.3 (14.7) 


16.4 (15.8) 


11 


12 


1.4 (0.5) 


16.0 (15.1) 


17.1 (16.2) 


16 


11 


2.0 (0.4) 


16.9 (15.3) 


18.0 (16.4) 



Scutum arm, and should be close to the corotation radius 
on the other side of the Galaxy. We have included the ex- 
pected extinction in the Northern and Southern Galactic 
hemispheres separately because they differ by a large mar- 
gin due to the presence of the Rho Ophiuchi cloud complex 
in the Northern hemisphere. Also included in Table [3] are 
the apparent magnitudes of B3 and B7 (reddened) stars at 
the indicated distances. It should be noted, however, that 
these reddening estimates are actually overestimates, since 
Gaia G is redder than Johnson- V making it less sensitive to 
interstellar absorption. 

For a B star with V < 15, we can infer from Table[3]that 
we will be able to observe stars with spectral type earlier 
than B3 as far as 12 kpc in both hemispheres and as far as 
14 kpc in the Southern hemisphere. On the other hand, stars 
with later spectral types (but earlier than B7) will be observ- 
able only in the Southern hemisphere but as far as 12 kpc. 
However, it will be possible to use B3 stars (and earlier) to 
probe distance as far as 18 kpc in the Southern hemisphere 
and 14 kpc in the Northern hemisphere, although with a 
loss of accuracy. In the same manner, it will be possible to 
use B7 stars (and earlier) as far as 14 kpc in the Southern 
hemisphere and as far as 12 kpc in the Northern hemisphere, 
although with less accuracy. 

However, we should remember that we are limited to 
distances of ~ 5 kpc and ~ 9 kpc for B3 and B7 stars, re- 
spectively, in terms of accurate distances determined from 
Gaia parallaxes, as was previously explained. Nevertheless, 
as was also pointed out, it should be possible to obtain accu- 
rate distances from follow-up spectroscopic analysis but this 
is within easy reach of modern medium-large telescopes. It 
is also important to note that this does not invalidate that, 
at least for a fraction of the sample (i.e. stars with spec- 
tral types earlier than B3), it will be possible to use direct 
distance determinations from Gaia. 

Thus, from the expected sample of ~ 2100 runaway 
stars 3 kpc above the Galactic plane, 525 are expected to be 
within 45° of the direction of the Galactic centre, if they are 
uniformly distributed. Hence, given the previous considera- 
tions, we expect to have a sample of ~ 250 in the Southern 
hemisphere covering a distance up to 18 kpc and ~ 175 in 
the Northern hemisphere covering a distance up to 14 kpc. 
The total number of stars should then be close to 400. 
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5 CONCLUSIONS 

The determination of the birthplaces of early-type stars is 
a feasible method to trace the position of the spiral arms in 
different time instants. We were able to estimate the pattern 
speed and the (present) phase angle of the Galaxy's spiral 
structure from a sample of local stars and obtained values 
(ftp = 20.3 ± 0.5 kms" 1 kpc" 1 and <f> = -53.1° ± 3.1°) that 
are consistent with previous estimates. Moreover, although 
these estimates are sensitive to systematic effects, in partic- 
ular biases introduced by stars born outside the spiral arms, 
these effects are not crippling and introduce an extra error 
of ~ 10 per cent. 

From the sample of high Galactic latitude runaway 
stars (sample B) we also derived a consistent estimate of 
the parameters (Q p — 21.9 ± 8.6 kms _1 kpc _1 and <f> = 
—63.5° ± 22.5°). However, given the nature of this sample 
(larger observational errors, small numbers and small spread 
in ages) the solution is not as well defined as in the case of 
the one obtained from sample A and has thus a larger un- 
certainty. 

More importantly, we have shown that this method has 
the potential to be used with data obtained by Gaia to re- 
trieve the birthplaces of distant runaway stars and hence 
it will be possible to trace not only the closest spiral arms 
but also the more distant ones, in particular the portions 
behind the Galactic centre, which are usually obscured by 
interstellar reddening. We have estimated that the number 
of runaway stars that will be available to trace these dis- 
tant portions of the spiral arms, in particular the far sides 
of the Sagittarius-Carina and Perseus arms, will be ~ 250 
to ~ 400, whereas the number of runaway stars available 
to trace spiral arms in closer regions of the Galaxy will be 
> 1000. Thus, considering the availability of the high qual- 
ity astrometric data obtained by Gaia in the near future, it 
will be possible to apply the method to a large number of 
runaway stars with high accuracy. 
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